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Abstract 

We consider the transverse-momentum [qx) distribution of Higgs bosons produced at 
hadron colliders. We use a formalism that uniformly treats both the small-gr and 
large-gr regions in QCD perturbation theory. At small qt {qr ^ Mh, Mh being the 
mass of the Higgs boson), we implement an all-order resummation of logarithmically- 
enhanced contributions up to next-to-next-to-leading logarithmic accuracy. At large 
9t {,<iT ^ Mh), we use fixed-order perturbation theory up to next-to- leading order. The 
resummed and fixed-order approaches are consistently matched by avoiding double- 
counting in the intermediate-g^ region. In this region, the introduction of unjustified 
higher-order terms is avoided by imposing unitarity constraints, so that the integral of 
the qT spectrum exactly reproduces the perturbative result for the total cross section 
up to next-to-next-to-leading order. Numerical results at the LHC are presented. 
These show that the main features of the distribution are quite stable with respect 
to perturbative QCD uncertainties. 
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The Standard Model (SM) of electroweak interactions has been spectacularly confirmed by 
experimental data. However the mechanism of mass generation remains to be understood. In 
its minimal version, the model predicts the existence of a scalar particle, the Higgs boson p], 
as a vehicle of electroweak symmetry breaking, but this particle has so far eluded experimental 
discovery. The LEP collaborations have put a lower limit on the mass Mh of the SM Higgs boson 
at about 114 GeV [2j, whereas fits to electroweak data prefer Mh<200 GeV at 95% CL [Hj. 
The next search for Higgs boson(s) will be carried out at hadron colliders, namely the Fermilab 
Tevatron g] and the CERN LHC 0. 

The main SM Higgs production mechanism at hadron colliders is the gluon fusion process. At 
leading order (LO), 0(ag), in the QCD coupling as this process occurs through a heavy-quark 
(top-quark) loop and, being a gluon-initiated process, it is expected to receive large radiative 
corrections. It is thus important to perform an accurate evaluation of higher-order QCD contri- 
butions, together with a reliable estimate of the associated theoretical uncertainty. 

The next-to-leading order (NLO) perturbative corrections to the total cross section for Higgs 
boson production via gluon fusion were computed in Refs. |U] (in the limit of an infinitely-heavy 
top quark) and [7| (including the dependence on the finite mass Mt of the top quark) and were 
found to be large (of the order of 80-100%), thus casting doubts upon the reliability of the 
perturbative expansion. In the last two years much effort has been devoted to improving the 
accuracy of the perturbative calculation. In the large- limit, the next-to- next-to-leading order 
(NNLO) contribution has been computed in Ref. |H] and still higher-order contributions have 
been evaluated in Ref. |H] by implementing soft-gluon resummation. Since these beyond-NLO 
corrections are moderate, the perturbative QCD predictions for the total cross section are under 
good control now. 

In this letter we consider a less inclusive observable, the transverse-momentum {qx) distribution 
of the Higgs boson. An accurate theoretical prediction of this observable at the LHC [S] can be 
important to enhance the statistical significance of the signal over the background and to improve 
strategies for the extraction of the Higgs boson signal. 

When studying the qx distribution of the Higgs boson in QCD perturbation theory, it is con- 
venient to start by considering separately the large-g^ and small-g^ regions. Roughly speaking, 
the large-gT region is identified by the condition qx^Mu- In this region, the perturbative series 
is controlled by a small expansion parameter, aQ,{M'jj)^ and calculations based on the truncation 
of the series at a fixed-order in as are theoretically justified^ In the small-gr region (g^r ^h)-, 
where the bulk of events is produced, the convergence of the fixed-order expansion is spoiled, since 
the coefficients of the perturbative series in as(M^) are enhanced by powers of large logarithmic 
terms, ln™'(M^/g|n). To obtain reliable perturbative predictions, these terms have to be system- 
atically resummed to all orders in as The fixed-order and resummed approaches have then 
to be consistently matched at intermediate values of gr, so as to avoid the introduction of ad-hoc 
boundaries between the large-g^ and small-g^ regions. 

Higgs boson production at large qx has to be accompanied by the radiation of at least one 
recoiling parton, so the LO term for this observable is of 0{a\). The LO calculation was reported 
in Ref. ^2]; it shows that the large- approximation works well as long as both Mh and qx are 

^We are not considering the extreme limit ^ where a resummation of enhanced perturbative terms is 
required |1U| . 
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smaller than Mj. Similar results on the validity of the large- approximation were obtained in the 
case of the associated production of a Higgs boson plus 2 jets (2 recoiling partons at large transverse 
momenta) In the framework of the large- Mf approximation, the NLO QCD corrections to 

the transverse-momentum distribution of the Higgs boson were computed first numerically ^1] 
and later analytically ^1^]. In the large-gr region, the overall effect of the NLO corrections to 
the qt distribution is of the same size as that of the NLO corrections to the total cross section. 

The method to systematically perform all-order resummation of logarithmically-enhanced 
terms at small qt is known ^2 IT7H^ (see also the list of references in Sect. 5 of Ref. |22j). 
To correctly take into account the kinematics constraint of transverse-momentum conservation, 
the resummation procedure has to be carried out in b space, where the impact parameter b is 
the variable conjugate to qt through a Fourier transformation. In the case of the Higgs boson, 
6-space resummation has been explicitly worked out at leading logarithmic (LL), next-to-leading 
logarithmic (NLL) [221 IM] and next-to-next-to-leading logarithmic (NNLL) j2S] level. The qt 
distribution is then obtained by performing the inverse Fourier (Bessel) transformation with re- 
spect to b. Various implementation formalisms [22 l26H3T| have been proposed to transform the 
resummed expressions back to qt space and to perform the matching with the fixed-order results 
at large qt- Phenomenological applications to the Higgs boson qt distribution have been presented 
in Refs. 1^ inH^I37j . by combining resummed and fixed-order perturbation theory at different 
levels of theoretical accuracy. 

In the following we use the formalism described in Ref. [HJ to compute the Higgs boson qt 
distribution at the LHC In particular, we combine the most advanced perturbative information 
that is available at present: NNLL resummation at small qt and NLO calculations at large qt- 
More details will be given elsewhere. 

We consider the collision of two hadrons hi and /i2 with centre-of-mass energy ^/s. According 
to the QCD factorization theorem (see Ref. [211 and references therein), the transverse- momentum 
differential cross section for the production of the SM Higgs boson can be written as 

—^igr, Mh,s) = dxi / dx2fa/hAxi,^l)fb/h2{x2,^l)-rYi9T,MH,s;as{iJ,R),fJ,%fJ,l), 

UQt ah ''^ '^'^ 

(1) 

where fa/h{x, fip) (a = q, q, g) are the parton densities of the colliding hadrons at the factorization 
scale ixp, ddab/dq^ are the partonic cross sections, s = X1X2S is the partonic centre-of-mass energy, 
and fiR is the renormalization scale. Throughout the paper we use parton densities as defined in 
the MS factorization scheme, and q;s(q'^) is the QCD running coupling in the MS renormalization 
scheme. 

The partonic cross section is computable in QCD perturbation theory and, as discussed above, 
it is evaluated by introducing the decomposition 



(res.) rl;ir^^^-^ 

+ (2) 



d^ab ^ da'^^ ■' ^ da^ 
dq^ dq^ dq^ 



The first term on the right-hand side contains all the logarithmically-enhanced contributions, 
ag/q^ ln"*Q^/gy, at small qx, and has to be evaluated by resumming them to all orders in as. 
The second term is free of such contributions, and can be computed by fixed-order truncation of 
the perturbative series. 
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The resummed component d&ac^'^ of the partonic cross section is written as 
— ^(gT,MH,s;as(Ai^),/i^,/i|) = - J db b Joibqr) Wadb, Mh, s; asifiu), IJ,%, IJ-I) , (3) 



dqj 

where Jo{x) is the Oth-order Bessel function. The factor W embodies the all-order dependence on 
the large logarithms L = lnM^6^ at large b, which corresponds to the gr-space terms In M^/g|. 
that are logarithmically enhanced at small qr (the limit gr Mh corresponds to Mnb ^ 1, 
because b is the variable conjugate to g^)- Resummation of these large logarithms is better 
expressed by defining the A^-moments Wat of W with respect to 2; = Mfj/s at fixed Mh'- 



Wac,N{b,MH;as{fi%),fi%,fil)= f dz z""-^ Wac{b, Mh, s = M^^ 

Jo 



/2;;as(/i^),/i^,/i|) . (4) 



The resummation structure of Wac, n can indeed be organized in exponential form as follows: 

Wjv(&, Mh; as{fil),fil, fil) = Hn [asif^l); Mfj/fil, Mfj/fil) 

X expig^iasifil), bMn; M^/fil, Mi//i|)} , (5) 

where the subscripts denoting the flavour indices are understood"'". 

All the large logarithmic terms ctgL"^ = In"* Mnb with 1 < m < 2n are included (actually, 
the complete dependence on b is included) in the form factor exp{^}. More importantly, all the 
logarithmic contributions to Q with n + 2 < m < 2n are vanishing. Thus, the exponent Q can 
systematically be expanded as 

gN{as,bMH;M'H/fi%,Ml/fil) = L g^^\asL) + {a^L; M^ / ^i^) 

+ asgPiasL;Ml/fil,Ml/fil) + ... , (6) 

where as = «s(Ai/j) and the functions g^^\asL) are defined such that g^'^^ = when asL = 0. 
Thus the term Lg^^^ collects the LL contributions a^L"''^^; the function g^"^^ resums the NLL 
contributions OgL"; g^^^ controls the NNLL terms OgL""^, and so forth. Note that in the expansion 
(jni) the logarithmic variable L has been replaced by 

L = \n{M'^bybl + l) , (7) 

where 60 = 2e~'^. In the resummation region Mnb ^ 1, the replacement is fully legitimate since 
L L. The reason for using L rather than L is discussed below. 

The function 'm. Eq. (0) does not depend on b and, hence, its evaluation does not require 
resummation of large logarithmic terms. It can be expanded in powers of as = '^s(/^/;) as 

Hjv(«s; Ml/iil, Ml/iil) = do «i [1 + g n^^\Ml/ii% Mil III) 

2 



+ ^%\Ml/^%Ml/^,l) + ...] , (8) 



■I" More precisely, "we are presenting the resummation formulae in a simplified form "which is valid "when there is 
a single species of partons. In general, the exponential is replaced by an exponential matrix "with respect to the 
flavour indices of the partons. 
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where cto = GF/(2887r-\/2) is the Bom level cross section in the large- Mj approximation, and 
Gf = 1.16639 X 10"^ GeV"^ is the Fermi constant. 

The 'finite' component da^^^'^ of the partonic cross section does not require resummation of 
large logarithmic terms either. We compute it as follows: 



d^ab 



(fill-) r^^ -, r^;v('"«'S.). 

(9) 



^2 



dqj^ L dqf \ f.o. L dq: 



d^ab 



f.o. 



The first term on the right-hand side is the usual perturbative series for the partonic cross section 
truncated at a given fixed order in as. The second term is obtained by truncating the resummed 
component in Eq. at the same fixed order in ag. The (small-g^) resummed and (large-gy) fixed- 
order approaches are thus consistently matched by avoiding double-counting in the intermediate-gT 
region. This procedure guarantees that the right-hand side of Eq. Q contains the full information 
of the perturbative calculation up to the fixed order specified by Eq. (jS)) plus resummation of 
logarithmically-enhanced contributions from higher orders. 

A few distinctive features of the formalism described so far require some comments. 

We implement perturbative QCD resummation at the level of the partonic cross section. In 
the factorization formula (Q), the parton densities are thus evaluated at the factorization scale 
\xfi as in the customary perturbative calculations at large qT- The central value of and /i^ 
has to be set equal to Mu-, the typical hard scale of the process, and the theoretical accuracy of 
the resummed calculation can be investigated as in fixed-order calculations, by varying \ip and 
[i^ around this central value. 

The variables L and L are equivalent to organize the resummation formalism in the region 
Mfffe ^ 1. The use of the variable L is inspired by the procedure introduced in Ref. [221 to 
deal with kinematical constraints when performing soft-gluon resummation in e^e~ event shapes. 
When Miib -C 1, L — » and exp{^} 1. Therefore, using the definition in Eq. ((7j), we avoid 
the introduction of all-order contributions in the small-6 region, where the use of the large-6 
resummation formalism is not justified. In particular, exp{^} = 1 at 6 = 0. This implies that 
the integral over q^ of da/dqx exactly reproduces the fixed-order calculation of the total cross 
section. Note that the bulk of the gr distribution is in the region qp^Mn- Since resummed and 
fixed-order perturbation theory controls the small-g^ and large-gr regions respectively, the total 
cross section constraint mainly acts on the size of the higher-order contributions introduced in the 
intermediate-gr region by the matching procedure. 

It is known pUH2Hl EUl 1^ that non-perturbative effects have an increasing role in the qx 
distribution as qx decreases. However, we do not include non-perturbative contributions. The 
main goal of the quantitative study presented below is to investigate the predictivity of QCD 
within a purely perturbative framework. In particular, we want to examine how the Higgs boson 
qT distribution is affected by perturbative QCD uncertainties, such as its dependence on scale 
variations and on higher-order contributions. 

The functions g]^\asL) and the coefficients H^^^ in Eqs. (jH)) and (jS)) can be expressed (see for 
instance Ref. [211) terms of perturbative coefficients known as A^"\ B^^\ C'-"^ [25 and H^'^^ j31j . 
In particular, g^'^^ depends on A^^'^\ g'^^ also depends on 5^^^ and A*^^) [22], g^N ^^^^ depends on 
H^^\ C(i) [211, ^^^^ PHITT)| and 7^!^^ depends on H^^^ and C^^\ u'ff also depends on H^^^ and 
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C^^). We also observe that the functions and receive additional contributions respectively 
from the LO and NLO anomalous dimensions that control the evolution of the parton densities. 
The NNLL coefficient A^^^ is not yet known. In the following we assume that its value is the same 
as the one [32] that appears in resummed calculations of soft-gluon contributions near threshold. 

(2) 

The coefficient 7i)y is not known in analytic form either. However, within our formalism we 
can exploit the property that the integral of the qt distribution exactly matches the fixed-order 
calculation of the total cross section. From the known NNLO result for the total cross section [Hj, 

(2,') iTT^n 

we thus extract Ti.)^ in (approximate) numerical form. As pointed out in Ref. [SJ, the coefficients 
B^"'\ C^"^ and H^"'^ cannot separately be defined without fixing a resummation scheme. Note, 
however, that the dependence on the choice of the resummation scheme cancels by recasting the 
resummed formulae in the form of Eq. (jSj): the functions g^^\asL) and the coefficients Tij^^in 
Eqs. dnj and (jH)) are explicitly resummation-scheme independent. 

The functions g^^\asL) are singular when A = (3oasL — > 1 (/5o is the first coefficient of the 
QCD /3-function). The singular behaviour is related to the presence of the Landau pole in the 
perturbative running of the QCD coupling as(g^). To properly define the b integration in Eq. Q, 
a prescription to deal with these singularities has to be introduced. Here we follow Ref. 021 and 
deform the integration contour in the complex b space, as an extension of the minimal prescription 
of Ref. m. 

In the following we present quantitative results at NLL+LO and NNLL+NLO accuracy. We 
implement Eqs. Q and 0. At NLL+LO accuracy, we compute da^"^^^'^ at NLL accuracy (we 
include the coefficient Ti^'* and the functions g]^^ and g^^^), and we match it with [(icrjf.o. evaluated 
at LO (i.e. at 0{a^)). At NNLL+NLO accuracy, we also include Ti^^^ and g^^^ in the resummed 
component and we evaluate [rfcrjf.o. at NLO (i.e. at 0{a^)). As for the evaluation of [(icr]f.o., 
we use the Monte Carlo program of Ref. ^l]- The numerical results are obtained by using the 
MRST2001 set of parton distributions |3S] and choosing Mh = 125 GeV. At NLL+LO we use LO 
parton densities and 1-loop as, whereas at NNLL+NLO we use NLO parton densities and 2-loop 
as- 

The NLL+LO results at the LHC are shown in Fig. ^ In the left-hand side, the full NLL+LO 
result (solid line) is compared with the LO one (dashed line) at the default scales = /i^ = Mh- 
We see that the LO calculation diverges to +oo as qx 0. The effect of the resummation is 
relevant below qt ~ 100 GeV. In the right-hand side we show the NLL+LO band that is obtained 
by varying fip = fJ'R between 1/2Mh and 2Mh- The scale dependence increases from about ±10% 
at the peak to about ±20% at qr = 100 GeV. The integral of the resummed curve is in good 
agreement with the value of the NLO total cross section evaluated with LO parton densities and 
1-loop as, the small difference being due to the (improvable) numerical precision of our code. 

The NNLL+NLO results at the LHC are shown in Fig. |21 In the left-hand side, the full result 
(solid fine) is compared with the NLO one (dashed line) at the default scales fip = f^R = Mh- 
The NLO result diverges to — oo as gr ^ and, at small values of gr, it has an unphysical peak 
(the top of the peak is above the vertical scale of the plot) which is produced by the numerical 
compensation of negative leading logarithmic and positive subleading logarithmic contributions. 
It is interesting to compare the LO and NLL+LO curves in Fig.^and the NLO curve in Fig. |21 At 
qt ~ 50 GeV, the qt distribution sizeably increases when going from LO to NLO and from NLO 
to NLL+LO. This implies that in the intermediate-gr region there are important contributions 
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Figure 1: LHC results at NLL+LO accuracy. 



that have to be resummed to all orders rather than simply evaluated at the next perturbative 
order. The qt distribution is (moderately) harder at NNLL+NLO than at NLL+LO accuracy. 
The height of the NNLL peak is a bit lower than the NLL one. This is mainly due to the fact that 
the total NNLO cross section (computed with NLO parton densities and 2-loop as), which fixes 
the value of the qt integral of our resummed result, is slightly smaller than the NLO one, whereas 
the high-gT tail is higher at NNLL order, thus leading to a reduction of the cross section at small 
qt- We find that the contribution of A^'^^ (recall that we are using an educated guess on the value 
of the coefficient A'^'^^) can safely be neglected. The coefficient IH!^ contributes significantly, and 
enhances the qt distribution by roughly 20% in the region of intermediate and small values of 
qt- The resummation effect starts to be visible below qt ~ 100 GeV, and it increases the NLO 
result by about 40% at qx = 50 GeV. The right-hand side of Fig. |21 shows the scale dependence 
computed as in Fig. [T] The scale dependence is now about ±6% at the peak and increases to ±20% 
at qt = 100 GeV. Comparing Figs. 1 and 2, we see that the NNLL+NLO band is smaller than the 
NLL+LO one and overlaps with the latter at qt ^ 100 GeV. This suggests a good convergence of 
the resummed perturbative expansion. 

We have considered perturbative QCD predictions for the Higgs boson qt distribution at 
the LHC. We have shown that the main features of the qt distribution are quite stable with 
respect to perturbative uncertainties (scale variations, inclusion of higher orders in the resummed 
expansion). More details about the formalism and our numerical results will be presented in a 
future publication, where we shall also consider the inclusion of non-perturbative contributions. 
Available studies |H3HTfj of non-perturbative contributions at the LHC estimate effects (at most) 
of the order of a few per cent when qt ^ 10 GeV. These effects are smaller than the resummation 
effects examined here. 

Acknowledgements 

We would like to thank Werner Vogelsang for discussions. S.C. wishes to thank Luca Trentadue 
for their early collaboration on Higgs boson production in hadronic collisions. 



6 




Figure 2: LHC results at NNLL+NLO accuracy. 
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